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Abstract. We present the results of quantitative spectral analyses of ten apparently normal B-type stars. These 
stars were found to be young massive B-type stars at distances of z=2.6 to 7.6 kpc from the galactic plane based 
on their positions in the (Teff, logg) diagram, normal abundance patterns and/or large projected rotational 
velocities. We discuss formation scenarios (runaway star scenarios versus a scenario for star formation in the 
halo) by comparing times-of-fiight and evolutionary time scales. For all stars (except SB 357 and HS 1914-1-7139) 
both the scales are similar indicating that the stars could have formed in the galactic disk and been ejected from 
there soon after their birth. Derived ejection velocities range from 130 km s~^ to 440 km s~^ and may be used to 
constrain models for ejection mechanisms. Using new proper motion measurements we show that PHL 346, which 
was considered the most likely candidate for a young B-type star born in the halo, can be explained as a runaway 
star from the galactic plane. 
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1. Introduction 

Main sequence B-type stars located far away from the 
galactic plane are a rare, albeit known phenomenon. In 
their pioneering paper, Greenstein & Sargent (1974) stud- 
ied faint blue stars at high galactic latitudes and classi- 
fied 25% of them as apparently normal OB-type stars at 
distances from the galactic plane of z=l-3kpc. More de- 
tailed studies have shown that many, but not all, of the 
apparently normal stars were in fact highly evolved low- 
mass stars. Surveys for UV excess objects (e.g., Palomar 
Green, Hamburg-Schmidt, Edinburgh Cape) have found 
many new candidates (e.g., Saffer et al 119971 RoUeston 



et al. 1999, Magee et al. 1998) and some may be known 
even in other galaxies (M 31, Smoker et al. 2000). Their 



properties, possible evolutionary histories and formation 
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mechanisms were reviewed by, e.g.,Tobin (1987), Keenan 
( |1992| ), an d Heb er et al. ( |1997| ). 

Tobin ( 1987 ) also discusses the problem that some 
highly evolved stars spectroscopically mimic massive 
stars almost perfectly. The most striking example is 
PG 0832-1-676 which has been analysed several times. 
Its abundance pattern is close to normal. Only re- 
cently, Hambly et al. (1996) were able to firmly establish 
slight underabundances and a very low projected rota- 
tion velocity. Combining both results they concluded that 
PG 0832+676 in fact is a highly evolved star. Abundance 
analyses as well as determinations of rotational veloci- 
ties are thus of essential importance for the verification 
of massive B-type star candidates. A high rotational ve- 
locity generally excludes a late evolutionary status of the 
star, as old, low-mass stars cannot rotate as fast as mas- 
sive stars. This fact was used, e.g., by Heber et al. (1995, 
HS 1914-^7139) and Schmidt et al. (|l99|, PG 0009+036) 
to identify massive B-type stars far from the Galactic 
plane from medium-resolution spectra. 

The massive B-type stars in galactic halos can be sep- 
arated kinematically into two different categories: those 
stars with a main sequence lifetime larger than the time 
they would need to travel from the plane to their present 
position and those with a main sequence lifetime too small 
to reach their current position assuming an acceptable ve- 
locity vertical to the galactic disk. The former ones are as- 
sumed to be born in the disk and thereafter ejected from it 
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(runaway stars), while the latter are supposed to be born 
in the halo (see Conlon et al. 198S, 199C, Hambly et al. 



1993 for more details). Since many years runaway stars 
are known to exist, whereas formation of massive stars 
in the halo has not yet been confirmed convincingly. The 
best studied candidate is PHL 346 (Ryans et al. 1996), a 
j3 Cephei star in the halo (Dufton et al. 19981 ). A search 
for coeval stars around PHL 346 (Hambly et al. 1996) 
met with limited success, since only one out of 16 A- and 
B-type stars around PHL 346 was found to have the ap- 
propriate spectral type and radial velocity. 

Calculations of galactic orbits are thus very important 
to determine the true nature of the stars, since they also 
allow to determine ejection velocities from the galactic 
disk. However, accurate proper motions are a prerequisite 
for such an analy sis. A big step forward has been achieved 
by ThejU et al. ( 1997 ) and the H ippar cos /Tycho mission 
(Ferryman et al. |1997| ; H0g et al. ^OOOl) . 

In this paper we present the analysis of new high- 
resolution spectra for 10 apparently normal B-type stars. 
Half of the sample are new discoveries whereas the other 
half has already been studied previously (Conlon et al. 
1992|, Ry ans et al. |1994 RoUeston et al. |1999| , Heber et 



al. 1995). For the latter the times-of-flight quoted in the 



literature appear to be larger than the evolutionary times 
indicating they might have formed in the halo. However, 
proper motions were not available rendering these esti- 
mates of the times-of-flight uncertain. Since proper mo- 
tion measurements became available recently for three of 
these stars, it was deemed necessary to reanalyse them 
from new high resolution spectra. 

2. Observations and Data Reduction 

We have obtained high resolution spectra for all pro- 
gramme stars using the HIRES spectrograph at the Keck 
I telescope, the FEROS spectrograph at the ESO 1.5m 
telescope, the FOCES spectrograph at the DSAZ 2.2m 
telescope and the CASPEC spectrograph at the ESO 3.6m 
telescope (see Table |l|) . Since the normalization procedure 
for the Echelle spectra is cumbersome for the rather broad 
Balmer lines (see below) it was deemed necessary to secure 
low resolution spectra at least for some of the programme 
stars to obtain independent estimates of the atmospheric 
parameters, in particular of the gravities. Appropriate 
spectra were obtained at Calar Alto and at ESO. Details 
are given in Table |l|. For PC 1511-^367, PC 1533+467 
and PC 1610-1-239 we have new low resolution spectra 
presented in this paper. These spectra were taken with 
the TWIN spectrograph at the DSAZ 3.5m telescope and 
reduced as described in Edelmann et al. ( ^001 ). The low 
resolution spectrum of PG 0122-)-214 were reduced as de- 



scribed in Moehler et al. (1997). 

For the Keck-HIRES spectra observed in 1998 the blue 
cross disperser was used and the spectra cover the blue 
spectral range (3600 to 5130 A). The sta ndard data reduc- 
tion as described by Zuckerman & Reid ( 1998| ) resulted in 
spectral orders that have a somewhat wavy continuum. 



Due to the merging of the higher Balmer lines the re- 
moval of the remaining waviness is difficult. We used the 
spectrum of H1504+65 (a very hot pre-white dwarf de- 
void of hydrogen and helium, Werner 1991), which was 
observed in the same night, for rectification of our spec- 
tra. Its spectrum has only few weak lines of highly ionized 
metals in the blue (3600~4480A) where the strong Balmer 
lines are found in the B-type stars. Therefore we normal- 
ized individual spectral orders 1 to 20 (3600-4480A) of the 
B-type stars by dividing through the smoothed spectrum 
of H1504+65. The remaining orders were normalized by 
fitting the continuum with spline functions (interpolated 
for orders 26 and 27 which contain H/3). Judging from the 
match of line profiles in the overlapping parts of neigh- 
boring orders this procedure worked extremely well. In 
1996 the red cross disperser of the Keck-HIRES spectro- 
graph was used and therefore the blue part of the spectrum 
shortward of 4200A was not recorded. Since no merging 
Balmer lines are present in this spectral range the wavi- 
ness of the spectrum could be removed by fitting the con- 
tinuum with spline functions. For orders containing broad 
Balmer lines the fit functions were interpolated between 
neighbouring orders. 

The FOC ES sp ectra were reduced as described in 
Pfeiffer et al. ( |l998|) with a software package developed by 
the Munich Group. For the FER OS data the MIDAS re- 
duction pipeline (Frangois, 1999 ) was used. For CAS PEC 
data the procedure described by Heber et al. ( 1986 ) was 
applied. The CASPEC and FEROS spectra were normal- 
ized in a similar way as described for the Keck HIRES 
spectra of 1996. 

Due to an error in the wavelength calibration no reliable 
radial velocity could be measured for HS 1914-1-7139. 



3. Atmospheric Parameters and Projected 
Rotational Velocities 

To derive atmospheric parameters (effective temperature, 
surface gravity, and photospheric helium abundance) and 
projected rotational velocities all Balmer lines and the He i 
lines AA 4026 A, 4388 A, 4438 A, 4472 A, 4713 A, 4922 A, 
5016A, 5048 A, 5678 A, in the observed spectr a were fitted 
with synthetic line profiles calculated from model atmo- 
spheres. 

We computed LTE model atmospheres using the pro- 
gram of Heber et al. (2000), which calculates plane par- 
allel, chemically homogeneous and fully line blanketed 
models, using the opacity distribution functions for metal 
line blanketing by Kurucz ( 1979| , ATLAS6). From these 
model atmospheres synthetic spectra were calculated with 
Lemke's versio n[| of the LINFOR program (developed 
originally by Holweger, Steffen, and Steenbock at Kiel 
University). The spectra include the Balmer lines H„ to 
H22 and the He i lines listed above and the grid covers the 



^ For a description see | http://a400.sternwarte.uni- 
|erlangen.de/^ai26 /Unfit /linfor.htnil| 
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Name 


Observation Date & 


Telescope & 


Resolution 


Wavelength Range 




Time (UT) 


Instrument 


(FWHM A) 


(A) 






Echelle spectra: 








PG 0122+214 


Jul 20, 1998 14:48 


Keck HIRES 


0.09 


3600 


- 5130 


PG 1511+367 


Jul 20, 1999 20:30 


Calar Alto 2.2m FOCES 


0.15 


3890 


- 6995 


PG 1533+467 


Jul 19, 1999 21:00 


Calar Alto 2.2m FOCES 


0.15 


3890 


- 6995 


PG 1610+239 


Jul 24, 1996 - 


Keck HIRES 


0.09 


4265 


- 6720 


PG 2219+094 


Jul 20, 1999 03:00 


Calar Alto 2.2m FOCES 


0.15 


3890 


- 6995 


PHL 159 


Sep 12, 1998 22:00 


Calar Alto 2.2m FOCES 


0.15 


3870 


- 6830 


PHL 346 


Oct 19, 1986 22:50 


ESO 3.6m CASPEC 


0.20 


4070 


- 5130 


SB 357 


Oct 1984 


ESO 3.6m CASPEC 


0.20 


4060 


- 5090 


HS 1914+7139 


Jun 05, 1996 13:08 


Keck HIRES 


0.09 


4280 


- 6720 


BD-15°115 


Sep 08, 2000 07:49 


ESO 1.5m FEROS 


0.09 


3630 


- 8860 






low resolution spectra: 






PG 0122+214 


Aug 07, 1990 02:32 


Calar Alto TWIN 3.5m 


3.5 


3875 


- 5010 


PG 1511+367 


Jul 18, 1999 20:49 


Calar Alto TWIN 3.5m 


2.9 


3350 


- 7700 


PG 1533+467 


Jul 18, 1999 21:00 


Calar Alto TWIN 3.5m 


2.9 


3350 


- 7700 


PG 1610+239 


Apr 11, 2001 


Calar Alto TWIN 3.5m 


1.0 


4100 


- 4950 


PG 2219+094* 


Jun 10, 1987 08:35 


ESO MPIA 2.2m B&C 


2.5 


4035 


- 4900 



* Moehler et al. (1990) 



range 11000 K< Tes < 40 000 K, 3.5 < logg < 6.5 and 
-4.0 < log ^ < -0.5 at solar metallicity. 

The fit procedure is based on a te st using the rou- 
tines developed by Bergeron et al. ( 1992 ) and SafFer et al. 
(1994) and modified by Heber et al. (1997) to derive also 
the rotational velocity. The theoretical spectra are con- 
volved with the instrument profiles (Gaussian with the ap- 
propriate instrumental FWHM) and a rotational profile. 
The fit program then normalizes theoretical and observed 
spectra using the same continuum points. Example fits 
for a rapidly rotating and a slowly rotating programme 
star are shown in Fig.^ for hydrogen and helium lines, 
whilst Figjl compares the metal line spectra of slowly ro- 
tating (PHL 159, BD— 15°115) stars and a rapidly rotat- 
ing star (PG 1533+467). The fit procedure was executed 
for all high and low resolution spectra and the results are 
listed in Table Formal fitting errors are very small 
for the high resolution spectra, on average AToff = lOOK, 
Alogg = 0.02. Systematic errors (e.g. continuum place- 
ment, uncertainties in line broadening theory) are cer- 
tainly larger and therefore dominate the error budget. We 
estimate errors in effective temperatures conservatively as 
5% and adopted an error of ±0.1 dex for the gravities. 
The fitting errors for the Helium abundance are on aver- 
age Alog(n(He)/n(H)) — 0.05. Since sharp Helium lines as 
well as broad Helium lines are well reproduced (see Fig. 
|l|) systematic errors due to He i line broadening theory 
appear to be small and we adopted an error of ±0.1 dex 
in all cases. 

For rapidly rotating stars the minimum is too 
poorly defined to allow a reliable determination of the 
He abundance simultaneously. Therefore, in a first step 




ZAMS 



log T,„ 



^ For the fit of the low resolution spectra the He abundance 
was kept fixed at the solar value and the projected rotation 
velocity derived from the high resolution spectra were used. 



Fig. 3. Positions of the programme stars (filled circles) in 
a (Tcff , logg) diagram with evolutionary tracks calculated 
by Schaller et al.(1992) for determining the masses and 
evolution times. 



the helium abundance was kept fixed at —1.00 (i.e. so- 
lar) for the fit procedure. In a second iteration step the 
helium abundance was determined by fitting the helium 
lines while keeping the TeS and log g fixed at those values 
determined in the first iteration step. For all stars (except 
HS 1914+7139) Str0mgren photometry is available, which 
allowed an independent determination of t he effe ctive tem- 
perature. We used the pr ogram of Moon ( 1985 ) as modi- 
fied by Napiwotzki et al. ( 1993 ) to derive the effective tem- 
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"log (g) = 3.81 +/- 0.1 
N(He)/N(H) = -0.97 +/- 0.008 
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T^„ = 18500 K +/- 900 K 
log (g) = 4.09 +/- 0.1 
N(He)/N(H) = -1.00 +/- 0.000 
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Fig. 1. Fit examples for a slowly rotating star (BD— 15°115, left hand side) and a rapidly rotating star (PG 1533+467, 
right hand side) 



Table 2. Atmospheric parameters and rotational velocities for the progamme stars as derived from high and low 
resolution spectroscopic data and comparison of these data with the effective temperatures calculated from Str0mgren 
photometry. The rotational velocities derived from the high resolution data were used to fit the low resolution spectra. 
*: Toff, logg fixed and He fitted. **: Helium is fixed at solar abundance, see text. 
For the low resolution data He is fixed at solar abundance. 



Name 



PG 0122+214 
PG 1511+367 
PG 1533+467 
PG 1610+239 
PG 2219+094 

PHL 159 
PHL 346 
SB 357 



BD-15°115 
HS 1914+7139 



Tcft 
(K) 



log( 



High Resolution 

^"6 N{H) 



(km s 



18 300 
16 100 

18 500 
15 500 

19 500 

18 500 

20 700 

19 700 



20 100 
17 600 



3.86 
4.15 
4.09 
3.72 
3.58 

3.59 
3.58 
3.90 



3.81 
3.90 



-0.98 
-1.16 
-0.94* 
-0.84* 
-1.00* 

-0.84 
-1.00 
-1.00** 



-0.97 
-0.99 



117 
77 
215 
75 
225 

21 
45 
180 



35 
250 



Low Resolution 



Teft 

(K) 



log(- 



18 700 
15 600 

17 700 
15 400 

18 200 



3.90 
4.20 
3.93 
3.69 
3.52 



18 100 



3.60 



Photometry 
Teft E(b-y) 
(K) 



18 500 (1) 

15 900 (1) 

17 700 (1) 

18 600 (1) 

16 700 (2) 

19 500 (3) 

20 900 (4) 
22 300 (7) 
19 700 (5) 
19 700 (8) 
19 800 (9) 

19 800 (5) 

20 200 (6) 



References: (1 ) Wesemael et al. \l99^; (2) Moon ey et al. ( ^00(\ ); (3) Moehl er et a l. il99(\ ); 
(4) Brown et al \l97^): ( 5) Hauck & Mermi lliod \l99^ ); (6) Kilkenny et al. \l97.^ ; 
(7) Kilkenny et al. \197\ ); (8) Kilkenny D. il99d^ ; (9) Graham et al. \l97^ ) 



0.0 
0.0 
0.020 
0.082 
0.037 
0.081 
0.025 
0.037 
0.052 
0.061 
0.037 
0.0 
0.0 
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Fig. 2. Wavelength range with strong N il and O ii lines to show examples for spectra of slowly rotating (bottom, 
middle) and rapidly rotating (top) stars. 



perature and the reddening and compare the photometric 
temperatures to the spectroscopic ones in Table ||. There is 
a good agreement between results from low and high reso- 
lution spectra and photometry, except for PG 1610+239, 
PHL 159 and PHL 346. The spectrum of SB 357 shows 
the presence of emission in H^ and H-y but not in Us. 
Therefore the effective temperature were obtained from 
Str0mgren photometry and the surface gravity from fit- 
ting the far wings of the hydrogen lines. The helium lines 
of this object were difficult to fit, but the observation is 
compatible with normal abundance and there is no indi- 
cation of emission in any of the helium lines observed. The 
parameters used for further analyses were taken from the 
high resolution spectra, because of the larger wavelength 
coverage and the excellent quality of the fits. In the case 
of PG 1533+467, however, the wavelength coverage of the 
low resolution spectrum is larger than that of the high 
resolution one and therefore we used the average. The fi- 
nally adopted parameters are listed in Table |[ Results are 
shown in a (Tcff, logg) diagram (Fig. ||). 



4. Chemical Abundances 

Seven programme stars (PG 0122+214, PG 1511+367, 
PG 1533+467, PG 1610+239, PG 2219+094, 



HS 1914+7139 and SB 357) display highly broad- 
ened lines (due to rotation, see Table ||). Only the 
strongest metal lines (e.g. C ii 4267A, Mg ii 448lA) could 
be identified. Therefore it was impossible to perform a 
detailed abundance analysis. 

The equivalent widths were measured employing 
the nonlinear least-squares Gaussian fitting routines in 
MIDAS with central wavelength, central intensity and full 
width at half maximum as adjustable parameters. For 
metal lines located in the wings of Balmer or helium lines 
an additional Lorentzian function is used to describe the 
line wings of the latter. 

Metal fines of the species C ii, C iii, N ii, O ii, Ne i, 
Mg II, Al II, Al III, Si II, Si III, P III, S ii, S iii, Ar ii 
and Fe ill were identified in the sharp-lined spectra of 
BD-15°115, PHL 159 and PHL 346. The atomic data for 
the analysis were taken from several tables: 



mm 



1. CNO from Wiese etal 

2. Fe from Kurucz ( |l992|) and Ekberg (|1993D 

3. Ne, Mg, Al, Si, S, P, Ar from Hirata et al. (|l995| ) 



The LTE abundances were derived by using the clas- 
sical curve-of-growth method and the LINFOR program 
(version of Lemke, see above). In this case the model 
atmospheres were generated for the appropriate values of 
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Fig. 4. LTE abundances (relative to l Her) and errors of the programme stars. Abundances derived from singly ionized 
elements are shown as filled circles and from doubly ionized ones as filled triangles. 
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0.0 



-0.5 



-1.0 



PHL 159 



C N Ne Mg Al Si P S Ar Fe 



Fig. 5. Like Fig. |4| LTE abundances (relative to l Her) 
for PHL 159. 

Abundances derived from neutral elements are shown as 
open circles, from singly ionized ones as filled circles and 
from doubly ionized ones as filled triangles. 



effective temperature, gravity and solar helium and metal 
abundance with the ATLAS9 program of Kurucz ( |1992| ). 



Then we calculated curves of growth for the observed 
metal lines, from which abundances were derived. Blends 
from different ions were omitted from the analysis. In the 
final step the abundances were determined from a de- 
tailed spectrum synthesis (using the LINFOR code de- 



scribed above) of all lines measured before. The results 
of the LTE abundance analysis and the r.m.s. errors for 
PHL 346 and BD-15°115 are shown in Table | and com- 
pared with other analyses and for PHL 159 in Table |[ 
Besides the statistical r.m.s. errors (given in Tables || and 
^) the uncertainties in Tea, logg and microturbulent ve- 
locity (see below) contribute to the error budget. In order 
to minimize the systematic errors we use the B-type star 
L Her as a comparison star. This star has been analysed 
by Hambly et al. ( |1997D . We redetermined the LTE abun- 
dances of L Her using the same atomic data, model atmo- 
sphere and spectrum synthesis code as for our programme 
stars and took the equivalent widths measured by Hambly 
et al. ( |1997| ). 

Our results for i Her agree to within 0.1 dcx with those 
of Hambly et al. ( |1997D except for C ii (0.12dex), Si ill 
(0.17 dex), S III (0.21 dex) and Fe ill (0.36 dex). In par- 
ticular our statistical error for Fe III is much lower than 



that of Hambly et al. (1997). These differences can be at- 
tributed to different oscillator strengths used. 
Results are given in Tables ^ and ^ and systematic errors 
are adopted for our programme stars as well. These errors 
are incorporated in the error bars plotted in Figs. ^ and ^ 
The determination of elemental abundances is interlocked 
with the microturbulent velocity ^. This can be derived if a 
sufficient number of lines of one ion can be measured over 
a wide range of line strengths. In our programme stars N 
II and O II lines are most suitable for this purpose since 
many lines of these ions can be identified. Microturbulent 
velocities of ^ = 8 km s~^ were found for PHL 159 and 
BD— 15°115, while a rather high value of ^ = 23 km s^^ 
was deduced for PHL 346. Our results for PHL 346 and 
BD— 15°115 are somewhat larger than those derived by 
Ryans et al. ( |1996D and Conlon et al. ( |1992D , see Table § 
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Table 3. LTE abundances of PHL 159 compared with l Her as a comparison star and the range of LTE abundances 
from 21 B-type stars analysed by Kilian (1994). For t Her the abundances are determined with the equivalent widths 
from Hambly et al. (1997). For b Her systematic errors due to uncertainties of atmospheric parameters have been 
determined in this work and are listed in parentheses. 



Element 


L Her 




PHL 159 
this paper 




B-type stars 
K94 


5 (km s-') 


5 




8 






He I 


10.78 




11.16±0.10 


(5) 




C II 


8.14 ±0.26 


(±0.21) 


8.17±0.35 


(5) 


8.02-8.95 


N II 


7.85 ±0.18 


(±0.14) 


7.85 ±0.36 


(9) 


7.48-8.30 


O II 


8.72 ±0.16 


(±0.19) 


8.98 ±0.17 


(17) 


8.24-8.65 


Ne I 






8.49 ±0.18 


(13) 




Ne II 










7.99-8.19'^' 


Mg II 


7.28 ±0.20 


(±0.12) 


6.67 


(1) 


7.02-7.68 


Al II 


6.18 










Al III 


6.31 ±0.12 


(±0.14) 


6.05 ±0.20 


(2) 


5.71-6.36 


Si II 


6.86 ±0.50 


(±0.18) 


6.38 


(1) 




Si III 


7.34 ±0.20 


(±0.21) 


7.64 ±0.24 


(4) 


6.73-7.65 


P III 


5.53 










S II 


6.99 ±0.19 


(±0.05) 


6.70 ±0.26 


(4) 




S III 


6.93 ±0.32 


(±0.17) 






6.23-7.48 


Ar II 


6.64 ±0.30 


(±0.08) 








Fe III 


7.33 ±0.09 


(±0.15) 


7.22 ±0.15 


(2) 


7.09-7.80 



References: K94=Kilian (1994) 

(1) = see text for a discussion of the Ne abundances 



Table 4. Comparison of LTE abundances for PHL 346 and BD— 15°115 with results from literature. The number of 
spectral lines used is given in brackets. Errors for the programme stars are statistical errors only. 



Element 


B-type stars 
K94 


PHL 346 
this paper 


PHL 346 
R96 


BD-15°115 
this paper 


BD-15°115 
C92 


C (km s"^) 




23 




15 




8 




5 




He I 




11.00±0.10 


(6) 


11.05 ±0.09 


(4) 


11.03±0.10 


(8) 


10.90 ±0.10 


(8) 


C II 


8.02-8.95 


8.14 ±0.32 


(2) 


8.16 ±0.29 


(3) 


8.11 ±0.20 


(4) 


8.00 ±0.30 


(3) 


C III 




8.33 


(1) 














N II 


7.48-8.30 


8.04 ±0.19 


(24) 


8.04 ±0.24 


(19) 


7.71 ±0.11 


(5) 


7.80 


(1) 


II 


8.24-8.65 


8.54 ±0.28 


(22) 


8.79 ±0.25 


(42) 


8.55 ±0.20 


(8) 


8.80 ±0.30 


(9) 


Mg II 


7.02-7.68 


7.14 


(1) 


7.37 


(1) 


7.13 


(1) 


7.10 


(1) 


Al II 
















6.30 ±0.10 


(2) 


Aim 


5.71-6.36 


6.27 ±0.07 


(3) 


6.20 ±0.06 


(3) 










Si II 




6.98 ±0.17 


(3) 


7.46 ±0.01 


(2) 










Si III 


6.73-7.65 


7.47 ±0.08 


(4) 


7.60 ±0.04 


(2) 


7.27±0.11 


(3) 


7.50 ±0.30 


(2) 


P III 




5.64 


(1) 


5.39 


(1) 










S II 




7.29 ±0.35 


(6) 






7.05 ±0.06 


(2) 


6.80 ±0.10 


(2) 


S III 


6.23-7.48 


7.02 


(1) 


7.20 ±0.09 


(7) 










Ar II 




6.72 ±0.09 


(2) 














Fe III 


7.09-7.80 


7.38 ±0.16 


(6) 


6.66 ±0.75 


(7) 


7.10 ±0.04 


(2) 


7.20 ±0.40 


(3) 



References: K94=KiUan (1994); R96^Ryans et al. (199t ); C92=Conlon et al. (199&) 



Remarkable is the large difference (« 1.0 dex) between 
the Si II and Si ill abundances. This has been found 
in several analyses of the comparison star l Her as well 
(Hambly et al. , 1997; 0.67 dex). In a differential analysis 



these systematic errors cancel to a large extent. NLTE ef- 



fects are small for all elements (<0.1 dex, Kilian 1994) ex- 



cept for Ne I. As demonstrated by Auer & Mihalas (1973) 
LTE calculations overestimate the neon abundance. They 
carried out NLTE calculations for Ne i in t Her and de- 



rived a neon abundance (close to solar) which is lower by 

0. 60 dex than our LTE result. Therefore our absolute Ne 
abundances are overestimated. The abundances of the pro- 
gramme stars with respect to t Her are plotted in Fig. ^. 

1. BD-15°115: 

All abundances are in good agreement with i Her (to 
within error limits) except for Mg and Fe which are 
underabundant by about 0.2 dex. Our results are in 
accordance with those of a previous analysis of the 
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star by Conlon ct al. (1992, see Table except for 
Si for which the authors derive a considerably larger 
abundance. 
PHL 346: 

Abundances of C, N, O, Al, Si, S, Ar and Fe arc in 
good agreement with those in t Her. P is enriched by 
about 0.15 dex, Mg is depleted by about 0.15 dex. Our 
results agree well with those derived by Ryans et al. 
( 1996D except for Si ii and Fe ill. For the latter our 
result has a much smaller error. 
PHL 159: 

Mg, Al and S are significantly depleted and O enriched 
by 0.3 dex, whereas the other elements are in good 
agreement with those of the comparison star. 



Spectral analyses of massive B-type stars in open clus- 
ters as well as in the field (e.g. Gies & Lambert 1992 



Table 5. Data of proper motion. The position angle is 
counted positive east of north. 



Name 


/i (mas/y) Position angle ° Reference 


PG 0122+214 
PG 1533+467 
PG 1610+239 
PG 2219+094 
BD-15°115 
PHL 346 


3.4 ± 2.3 234 ±71 1 
16.8 ± 3.6 326 ± 1 

8.1 ± 4.0 150 ±16 1 

6.2 ± 3.6 194 ±47 1 
9.0 ± 2.7 92 ± 10 2 
8.9 ± 3.1 144 ±20 3 



References: (1) Thefl l et a l. \199'^ 
(2) Perryman et al. (1991); 



(3) Tycho-2 catalog, H0g et al. (200C ) 



1.11 



a 



IQO.iVo [kpc] 



Kilian |T994[ Cunha & Lambert |l994| ) have revealed con- ^^lere is the stellar mass in Mg, g is the gravity in 



siderable variations of metal abundances from star to star ^ model atmosphere flux at the stellar sur 

face in units of 10^ erg cm^^s^^A^^ and Vq is the dered- 
dened apparent visual magnitude. 

6. Kinematics 

6.1. Radial Velocities and Proper Motions 

Radial velocities of the slowly rotating programme stars 
were derived from the lineshift of metal lines. For the 
rapidly rotating stars only the Balmer and He i lines could 
be used. Radial velocities obtained this way were then cor- 
rected to heliocentric values and the results are listed in 
Table | The error of the velocities estimated from the 
scatter of the velocities derived from individual lines is 
about 3-11 km s~^.Our measurements agree to within 
error limits with previous estimates (see Table ||). 
Proper motions were taken from literature and are listed 
in Table |. 

6.2. Times-of-Flight and Ejection Velocities 

The times-of-flight, which the stars need to reach their cur- 
rent halo positions from the galactic disk, were calculated 
with the program 0RBIT6 developed by Odenkirchen &; 
Brosche ( 1992| ). This numerical code calculates the orbit of 
a test body in the Galactic potential of Allen & Santillan 
(1991). The complete set of cylindrical coordinates is in- 
tegrated and positions and velocities are calculated in 
equidistant time steps. The input for this program ver- 
sion are equatorial coordinates, distance d from the sun, 
heliocentric radial velocities and observed absolute proper 
motions. Values for proper motions are given in Table ^ 
The proper motions for PHL 159, PG 1511+467, SB 357 
and HS 1914+7139 were set to zero, because no measure- 
ments are available. We followed the orbits backwards in 
time (time steps of 0.01 - 0.1 Myr). The time of passage 
through the galactic disk (= change of sign in z-position 
relative to the disk) defines the time-of-flight Tflight- The 
velocity at the time of first crossing of the galactic plane is 



(even within an open cluster). Kilian (1994) carried out 
spectral analyses of 21 B-type stars in two open clusters 
and in the field and determined abundances of C, N, O, 
Ne, Mg, Al, Si, S, and Fe. We compare our results for 
PHL 346, BD-15°115, PHL 159 and l Her to her LTE 
results in Tables || and ^. Since her programme stars are 
somewhat hotter than ours, the Ne abundance is based on 
Ne II lines, whereas we had to use Ne i lines. Correcting for 
the significant NLTE effect on Ne I (0.56 dex, see above) 
the neon abundance of PHL 159, the only programme star 
for which it has been measured, is found to be consistent 
with Kilian's distribution. The abundances we derived for 
aU metals of PHL 346, BD-15°115 and t Her lie well 
within Kilian's distribution indicating that they are bona 
fide main sequence B-type stars. For PHL 159, however the 
O II abundance is higher and the Mg ii abundance lower 
than in Kilian's distribution, whereas the other metals 
are consistent with that distribution. Therefore PHL 159 
might either be a massive B-type star with rather peculiar 
abundances of the elements O and Mg or an evolved, low 
mass B-type star that mimics a massive B-type star quite 
closely. 

5. Masses, Distances and Evolutionary Times 

The derived atmospheric parameters were compared to 
two sets of evolutionary tracks (from the Geneva group, 
Schaller et al. 1992, and the Padua group, Salasnich et al. 
200C) to estimate stellar masses and evolutionary times 
Tevoi (see Fig. ^) by interpolation. The results do not de- 
pend on the model grid used. Derived masses differ by less 
than 0.1 Mg and evolutionary lifetimes by less than 4 Myr 
(except PG 1610+239: 17 Myr). Errors for the evolution- 
ary lifetimes in Table |^ include errors propagated from 
uncertainties in atmospheric parameters as well as from 
the use of the two model sets. 

The distance has been calculated from mass, effective 
temperature, gravity and the dereddened apparent mag- 
nitude of the stars: 
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regarded as the ejection velocity Vc and is also calculated 
by the program 0RBIT6. 

Results for all parameters of the programme stars (ef- 
fective temperature, gravity, radial velocity, ejection veloc- 
ity, mass, distance, age and time-of-flight) are summarised 
in Table |^. For the origin of the stars (see next section) 
the ages (Tovoi) and the times-of-flight (Tfiight) are im- 
portant. We improved Tflight for BD-15°115 and PHL 346 
for which proper motion measurements have become avail- 
able recently. For BD— 15°115 we derive a somewhat lower 
Tfiight than Conlon et al. (1992) and find Tflight to be con- 
sistent with Tovoi to within our error limits. For PHL 346 
we confirm that Tflight is slightly larger than Tcvoi, but 
given the error limits this is insignificant. 
For PG 2219-1-094 we find Tflight and Tovoi to be lower 



than derived by RoUeston et al. ( |1999| ) and Tflight to be 
consistent with Tevoi- For PG 1610-^239 the time-of-flight 
is poorly constrained and only a lower limit could be de- 
termined which is consistent with the estimate of the evo- 
lutionary time. 

7. Discussion 

We have carried out quantitative spectral analyses of ten 
apparently normal B-type stars. Their positions in the 
(Tefi, logg) diagram are consistent with models for main 
sequence stars. SB 357 shows emission in and H^, 
which confirms its classification as a Be star (Heber & 
Langhans 1986, Kilkenny 1989). Seven of the stars have 
rotational velocities > 70 km s^^ making detailed abun- 
dance analyses impossible. Mostly normal abundances 
with respect to l Her were determined for BD— 15°115, 
PHL 159 and PHL 346. The Mg and O abundances of 
PHL 159 are significantly different from th e com parison 
star and other normal B-type stars (Kilian, 1994). 

Calculated orbits based on measurements of radial ve- 
locity and proper motion allowed to determine times-of- 
flight from the galactic plane to their present position. 

7.1. Runaway stars 

Times-of-flight for PG 1511+367 and PG 1533-K467 are 
smaller than the evolutionary times, indicating that these 
stars have been formed in the galactic plane and were then 
ejected (runaway stars). The times-of-flight are similar to 
the evolutionary times for PG 0122-f 214, PG 2219-f 094, 
PHL 159, BD-15°115, and PG 1610-1-239, which implies 
that the stars could also have formed in the galactic disk 
and were then ejected very soon after their birth. Ejection 
velocities for all programme stars range from 130km s~^ 
to 440 km s~^. 

Three mechanisms for the production of runaway stars 
have been proposed in the literature: 



i) In the binar y supernova scenario (Zwicky, 1957 ; 
Blaauw, 1961 ) the runaway star receives its velocity 
after a supernova explosion in a massive close binary. 
After the explosion the binary occasionally dissociates 



and the secondary is traveling with a velocity compa- 
rable to its pre-explosion orbital velocity. Calculations 
by Iben & Tutukov, (1996) indicate that the runaway 
star can gain a velocity of 100 to 200km s^^. 
ii) In the dynamical ejection scenario (Poveda et al 

Gies & Bolton, 1986) the runaway star gains its veloc- 



1967 



ity through a dynamical interaction with one or more 
stars. The most efficient encounter is that between 
two close binaries in a stellar association or open clus- 
ter, which in most cases results in the ejection of two 
runaway stars and one eccentric binary ( Hoffer 1983; 



Mikkola 1983 ). Calculations by Leonard (1991) show 



that even velocities in excess of 1000 km s^ can be 
gained in rare cases, 
iii) It has been conjectured that star formation can be 
triggered by the interaction of the gas of the galactic 
disk with an infalling satellite dwarf galaxy. The mo- 
mentum transferred to the gas results in a significant 
velocity component of the newly born stars perpen- 
dicular to the galactic plane. This scenario has been 
supported by the discovery of the encounter with the 
Sagittarius dwarf galaxy (Ibata et al. 1994). 



The availability of precise milliarcsecond astrometry 
for nearby stars through the Hipparcos satellite and pul- 
sar astrometry and timing measurements have recently 
demonstrated that both the binary supernova scenario 
and the dynamical ejection scenario are viable. By cal- 
culating orbits for runaway stars, pulsars and open clus- 
ters it recently became possible to associate runaway stars 
and pulsars with their nascent clusters (Hoogerwerf et al. , 
^OOOj 120011) . 

Our programme stars are too far away and their space 
motions are therefore not known accurately enough to al- 
low to identify their relation to a young cluster or associa- 
tion in the galactic plane. However, the ejection velocities 
determined for our programme stars may be important 
to identify the mechanism which led to their ejection from 
the galactic plane, once reliable theoretical predictions be- 
come available for the different scenarios discussed above. 
Six stars have escape velocities exceeding 300 km s~^ 
which seems too large to be achievable by the binary su- 
pernova scenario. 

7.2. Stars born in the galactic halo? 

Four stars in our sample have been proposed in the litera- 
ture as candidates for B-type stars formed in the galactic 
halo because their times-of-flight were found to be consid- 
erably larger than the evolutionary time scale s (Co nlon 



et al., 1992, Keenan et al., 1986, Heber et al., 1995). As 



discussed above our new analysis of BD— 15°115 demon- 
strates that its time-of flight is consistent with the evolu- 
tionary time. Hence it could be runaway star, too. 

For SB 357 and HS 1914-^7139 the times of flight are 
more than twice as large as the evolutionary times, which 
would make formation in the disk unlikely. However, their 
times-of-flight are uncertain due to the lack of proper mo- 
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Table 6. Physical and kinematic parameters of the programme stars. Five stars have been analysed previously. The 
results of these analyses from literature are given for comparison. The errors of the evolutionary times are formal 
errors introduced by the errors in Toff and logg. 



Name 


Tcff 
K 




Vrad 
km S"""^ 


Vc 

km 


M 

M0 


d 
kpc 


z 
kpc 


Tffight 
Myr 


Tcvol 

Myr 


FG 0122+214 


18 300 


3.86 


26+5 


290 


6.7 


9.6 


6.2 


51 + 24 


35 + 6 


PG 1511+367 


16100 


4.15 


102+11 


300: 


4.8 


3.8 


3.2 


24: 


34 + 7 


PG 1533+467 


18100 


4.00 


33+6 


440 


6.0 


3.0 


2.4 


20 + 4 


33 + 5 


PG 1610+239 


15 500 


3.72 


91+10 


130 


5.8 


8.4 


5.9 


> 62 


54 + 10 


PHL 159 


18 500 


3.59 


88+3 


320: 


8.0 


5.3 


3.2 


31: 


28 + 2 


PG 2219+094 


19 500 


3.58 


-24+9 


220 


8.7 


9.8 


6.1 


43 + 22 


27 + 2 


(1) 


17 900 


3.60 


-7 


- 


7.5 


- 


- 


67 


41 


BD-15°115 


20100 


3.81 


93+4 


410 


8.0 


4.9 


4.8 


30 + 5 


26 + 4 


(2) 


19 500 


3.50 


94 




10.0 






47 


20 


HS 1914+7139 


17 600 


3.90 




330: 


6.2 


14.9 


6.0 


91:* 


39 + 6 


(3) 


18 000 


3.75 


-39 




6.5 - 10.0 


16 - 18.4 








PHL 346 


20 700 


3.58 


63+4 


350 


9.9 


8.7 


7.4 


27 + 7 


19 + 2 


(4) 


22 600 


3.60 


66+10 




13.0 




8.7 
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SB 357 


19 700 


3.90 


58+10 


230: 


7.4 


7.9 


7.8 


61: 


26 + 4 


(2) 


19 000 


3.70 


54 




8 




9.0 


64 


25 



References: (1) Rolleston et al. jWO^ ; (2) Conlon et al. \l99i) : (3) Heber et al. !\199^ ; (4) Keenan et al. \l98d^ 
* based on the radial velocity derived by Heber et al. (1995) 
: uncertain due to the lack of proper motion data 



tion measurements. Such data are urgently needed before 
any firm conclusions can be drawn. Therefore we are re- 
luctant to regard these stars as born in the halo. 

PHL 346 has been proposed as a candidate massive 



B-type star born in the halo (Ryans et al. 1996, Hambly 
et al. 1996| ) . Based on the new Tycho proper motion mea- 
surement, our analysis indicates that Tf^ght is marginally 
larger than Tcvoi and PHL 346 can be a runaway star, too. 

Hence no conclusive candidate for a young massive B- 
type star formed in the halo remains in our sample. Proper 
motions for the four stars lacking any measurement should 
urgently be determined. 
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